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Abstract 

The dynamics of spinning particles in curved space-time is discussed, emphasizing the hamiltonian formu¬ 
lation. Different choices of hamiltonians allow for the description of different gravitating systems. We give 
full results for the simplest case with minimal hamiltonian, constructing constants of motion including 
spin. The analysis is illustrated by the example of motion in Schwarzschild space-time. We also discuss 
a non-minimal extension of the hamiltonian giving rise to a gravitational equivalent of the Stern-Gerlach 
force. We show that this extension respects a large class of known constants of motion for the minimal 
case. 
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1 Spinning-particle dynamics 

The dynamics of angular momentum and spin of gravitating compact bodies has been a 
subject of great interest and intense investigation since the early days of relativity theory 
[DElEllllElEUZUSlEUinUIIlE]; for recent overviews see mmm- As argued in [16] 
there are two complementary approaches to the subject. One approach starts from the 
covariant divergence-free energy-momentum tensor of matter, which makes it possible to 
keep track of aspects of the structure of the body. The energy-momentum vector and the 
angular-momentum tensor can be constructed by computing integrals of components of 
the energy-momentum tensor and their first moments over the volume of the body, using 
suitable boundary conditions. Equations of motion for these quantities are then derived 
by applying the conservation law for the energy-momentum tensor of matter 0E10J. 

The other approach is to construct effective equations of motion for point-like objects, 
which is an idealization of a compact body, at the price of neglecting details of the internal 
structure by assigning the point-like object an overall position, momentum and spin. This 
is also known as the spinning-particle approximation, and is used for the semi-classical 
description of elementary particles as well. A large variety of models for spinning particles 
is found in the literature PinillSlIISlEniEIlESlESlEaESlISSlEZl. 

In this letter we take the second point of view for the description of spinning test 
masses in curved space-time, using an effective hamiltonian formalism similar to the one 
introduced in ref. [28]. One of the advantages of this description is that it can be applied 
to compact bodies with different types of spin dynamics, such as different gravimagnetic 
ratios. In this way specific aspects of the structure can still be accounted for. 


2 Covariant phase-space structure 

Hamiltonian dynamical systems are specified by three sets of ingredients: the phase space, 
identifying the dynamical degrees of freedom, the Poisson-Dirac brackets defining a sym- 
plectic structure, and the hamiltonian generating the evolution of the system with given 
initial conditions by specifying a curve in the phase space passing through the initial point. 
The parametrization of phase-space is not unique, as is familiar from the Hamilton-Jacobi 
theory of dynamical systems. Changes in the parametrization can be compensated by 
redefining the brackets and the hamiltonian. A convenient starting point for models with 
gauge-field interactions is the use of covariant, i.e. kinetic, momenta rather than canoni¬ 
cal momenta; see [29] and references cited there for a general discussion, and [28] for the 
application to spinning particles. 

The spin degrees of freedom are described by an antisymmetric tensor E Mi/ , which can 
be decomposed into two space-like four-vectors by introducing a time-like unit vector u\ 
= — 1, and defining 
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By construction both four-vectors S and Z are space-like: 


A'X, = 0, = 0. 


( 2 ) 
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In the following we take u to be the proper four-velocity of the particle. Then S is the 
Pauli-Lubanski pseudo-vector, from which a magnetic dipole moment can be constructed, 
whilst the components of Z, which will be referred to as the Pirani vector, can be used 
to define an electric [3HJ or mass dipole moment |311 32], Observe that we can invert the 
relations ([TJ) to write 

YT = - ^e‘ WKX u K S x + vPZ v -v^Z*. (3) 

V~9 

Therefore, if the Pirani vector vanishes: Z — 0 [33], the full spin tensor can be reconstructed 
from S. However, in non-flat space-time this is generally not the case. 

The full set of phase-space co-ordinates of a spinning particle thus consists of the 
position co-ordinate x 0 the covariant momentum and the spin tensor XP", with anti¬ 
symmetric Dirac-Poisson brackets 

{x M , ir v } = 6£, {vr M , 7r„} = ^ E K ' A /?. reA/w , 

= < 4 ) 

+ g vX S'**. 

The brackets imply that 7r represents the generator of covariant translations, whilst the 
spin degrees of freedom E generate internal rotations and Lorentz transformations. It is 
straightforward to check that these brackets are closed in the sense that they satisfy the 
Jacobi identities for triple bracket expressions. Thus they define a consistent symplectic 
structure on the phase spac^} 

To get a well-defined dynamical system we need to complete the phase-space structure 
with a hamiltonian generating the proper-time evolution of the system. In principle a large 
variety of covariant expressions can be constructed; however if we impose the additional 
condition that the particle interacts only gravitationally and that in the limit of vanishing 
spin the motion reduces to geodesic motion, the variety is reduced to hamiltonians 

H = H 0 + H s , H 0 = ^~ tt^, (5) 

2m 

where H^, = 0 whenever E^ = 0. In this letter we focus first on the dynamics generated 
by the minimal hamiltonian Hq. However, we also consider an extension with mi 

Hz = ^ ( 6 ) 

The choice of hamiltonians can be enlarged further by including charges coupling the 
particle to vector fields like the electromagnetic held [El] [30]. 


1 We have not found this complete set of brackets in curved space-time in the literature. However, other 
sets of brackets have been proposed [3] based on a larger set of degrees of freedom, some of which are 
subsequently removed by supplementary constraints. 
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3 Equations of motion 

Eqs. Q and (J5| specify a complete and consistent dynamical scheme for spinning particles. 
Note that the choice of hamiltonian is fixed by further physical requirements, and can differ 
for different compact objects. In that sense the hamiltonian is an effective hamiltonian, 
suitable to describe the motion of various types of objects in so far as the role of other 
internal degrees of freedom can be restricted to their effects on overall position, linear 
momentum and spin. 

The simplest model is obtained by restricting the hamiltonian to the minimal geodesic 
term Ho- By itself this hamiltonian generates the following set of proper-time evolution 
equations: 

= {x^, H 0 } => = rng^x 1 ', (7) 

stating that the covariant momentum n is a tangent vector to the world line, proportional 
to the proper four-velocity u — x. Next 

Ti = {*», H 0 } => D t - x x Y x fn u = Z kX R kX " tt u , (8) 

which specifies how the world line curves in terms of the evolution of its tangent vector. 
Finally the rate of change of the spin tensor is 


t> w = {E^, Ho} =* D r YZ V = t' llJ + x x Y x fZ KV + x A T 


\k 


E^ K = 0. 


(9) 


In these equations the overdot denotes an ordinary derivative w.r.t. proper time r, whereas 
D t denotes the pull-back of the covariant derivative along the world line x^[r). By sub¬ 
stitution of eq. ([7]) into eq. ([8]) one finds that 


Dix >l = x^ + Yj l x x x v 


Xu 


= —s nX RJ l u x v 
2m kA v 


( 10 ) 


which reduces to the geodesic equation in the limit E = 0. The world line is the solution 
of the combined equations (10) and ([9]) satisfying some initial conditions. This world line 
is a curve in space-time along which the spin tensor is covariantly constant. It has been 
remarked by many authors that the spin-dependent force (|8| exerted by the 

space-time curvature on the particle is similar to the Lorentz force with spin replacing the 
electric charge and curvature replacing the electromagnetic held strength. In this analogy 
the covariant conservation of spin along the world line is the natural equivalent of the 
conservation of charge. 

Even though the spin tensor is covariantly constant, this does not hold for the Pauli- 
Lubanski and Pirani vectors S and Z individually. Indeed, due to the gravitational Lorentz 
force 

D T S» = 4 m \— g^ KA E KA E a PR aPvp u p , 

( 11 ) 

d t z^ = ^e^e^^x, 
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where is the linear expression in terms of and Z M given in eq. (J3]) . We observe that 
the rate of change of both spin vectors is of order C[E 2 ]. In particular, as Z is not conserved 
in non-flat space-times the condition Z — 0 cannot be imposed during the complete motion 
in general. Indeed, the evolution of the system is completely determined by eqs. ®|@. 
and leaves no room for additional constraints. 

We close this section by remarking that the gravitational Lorentz force for unit mass 
1/2 T, kX R k £' 1/ u 1 ' can be interpreted geometrically as the change in the unit vector tA gener¬ 
ated by transporting it around a closed loop with area projection in the aA-aWplane equal 
to E kA . 


4 Conservation laws 

By construction the time-independent hamiltonian represented by (J5]) , (Jb]) is a constant 
of motion for the spinning body, irrespective of the specific geometry of the space-time 
manifold. In particular for the minimal geodesic hamiltonian Hq we have 


° 2 


Another obvious constant of motion is the total spin: 


/ = 2 9**9 a,S kA E^ = + Z^. 


( 12 ) 


(13) 


In addition, there may exist conserved quantities E) resulting from symmetries of 

the background geometry, as implied by Noether’s theorem OHDIES!. They are solutions 
of the generic equation 


{J,H 0 } = -<rn v 

m 


— + r - 
dx » T ^ A 


- dJ | 1 ya gp dJ | r 9J 

K c>7r A + 2 S apXfl dTT X + ^ a <9E* A 


= 0. (14) 


It follows that any constants of motion linear in momentum [36] are of the form 

J = (15) 

with 

V ^ol v + = 0, V A /V = R^\CXk- (16) 

These equations imply that a is a Killing vector on the space-time, and f3 is its anti¬ 
symmetrized gradient: 

/V = ^ (V^cty - VjyCx^). (17) 

Similarly constants of motion quadratic in momentum [37] are of the form: 

J = \ + 1 (18) 
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where the coefficients have to satisfy the ordinary partial differential equations 

VA OtflV “I” V [lOtyX ~\~ V 0 , 

Vh/3k\ v + VufinXn = a vp + R K \u a ppi ( 19 ) 


V P 7m-A — R^p (3 k\ a + R K \p ft; 


liver- 


Thus a is a symmetric rank-two Killing tensor, and the coefficients (/3, 7 ) satisfy a hierarchy 
of inhomogeneous Killing-like equations determined by the ap V . In the case of Grassmann- 
valued spin tensors E^" = ii(j ffi/’ 1 ' the coefficient 7 is completely anti-symmetric and the 
equations are known to have a solution in terms of Killing-Yano tensors 


The constants of motion (15) linear in momentum are special in that they define a Lie 
algebra: if J and J’ are two such constants of motion, then their bracket is a constant of 
motion of the same type. This follows from the Jacobi identity 


{{J, J'} , H 0 } = {{J, H 0 } , J'} - {{J', H 0 } , J} = 0. (20) 

Thus, if {ei}} =1 is a complete basis for Killing vectors: 

effi = edef, eJV.e? - e?V„eJ = /■/<, 
the constants of motion define a representation of the same algebra: 

J i = e^ + ^V fi e iv ^ =► {•/,../,} ( 21 ) 

Evidently such constants of motion are helpful in the analysis of spinning particle dynamics 

[mum Ei- 

5 Schwarzschild space-time 

The dynamics of spinning bodies can be illustrated by the motion in a static and spherically 
symmetric Schwarzschild space-time, for which the hamiltonian Hq in Droste co-ordinates 
is given by 


2mHo = 


1 _ 2 M 


< + 1 


2 M 


"2i 1 2 2 , 2 - 2/12 

7T T V Tin T T Sill 0 IT 


( 22 ) 


The space-time manifold admits four Killing vectors, for time-translations and rotations. 
They give rise to the conservation of kinetic energy: 


— E 


v , M ytr 
n t ~\—v E , 


(23) 
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and angular momentum: 


J\ — — sin (pne — cotan 9 cos p tt v 

—r sin p E r ’ e — r sin 9 cos 9 cos p Ti rtp + r 2 sin 2 9 cos p E 61 ^, 


J 2 = cos p 7r e — cotan 9 sin p 

+r cos p Y7 e — r sin 9 cos 9 sin p + r 2 sin 2 9 sin p E 6 ^, 


( 24 ) 


J 3 = vry, + r sin 2 9 E r ^ + r 2 sin 9 cos 9 E 6 ^. 

It is straightforward to check that these satisfy the usual algebra of time-translations and 
spatial rotations: 

{E,Ji} = 0, {Ji, Jj} = EijkJk- (25) 

As usual, the conservation of total angular momentum and the spherical symmetry of the 
space-time geometry allow one to take the angular momentum J as the direction of the 
z-axis, such that 

J — (0, 0, J). (26) 

For spinless particles, for which the angular momentum is strictly orbital, this implies that 
the orbital motion is in a plane perpendicular to the angular momentum 3-vector; with 
our choice of the z-axis this is the equatorial plane 9 = tt /2. 

In the presence of spin the result no longer holds in general, as the precession of spin can 
be compensated by precession of the orbital angular momentum, resulting in a non-planar 
orbit |30|. However, one can ask under which conditions planar motion is still possible. As 
in that case the directions of orbital and spin angular momentum are separately preserved, 
it means that necessary conditions for motion in the equatorial plane are 


J\ — J2 — 0 , 7 Tg — 0 , 


(27) 


and therefore also 

E r0 = E^ = 0. (28) 

Furthermore the absence of acceleration perpendicular to the equatorial plane expressed 
by D t 7Tq = 0 implies that 

E ie = 0. (29) 


Thus planar motion requires alignment of the spin with the orbital angular momentum; it 
is straighforward to show that the reverse statement also holds [3H ! 12 . 

In terms of the four-velocity components we are now left with relevant constants of 
motion 


E = m 




u 


M tr 
r 2 ^ ’ 


(30) 


and 


J = mr 2 iE + rTi rip , 


(31) 


6 



in addition to the hamiltonian constraint 


1 - 


2 M 


u t2 = 1 + 


u 


r 2 


| 2 M 


i 2 w2 

T- r iv , 


(32) 


and the conservation of total spin /, or equivalently: 

1 T y tr 2 

1 7 ^ | 




r 2 -[_ _ 2M ' _ 2M y 


(33) 


2 Mj 2 ' 

These equations show, that once the orbital velocities are known, all the non-vanishing 
spin components can be calculated from eqs. (30), (31) and (33). 

The simplest type of planar orbit is the circular orbit r = R = constant, u r = 0. In 
this case the symmetry of the orbit implies that (td, vJ) are constant in time, and that 
T, tip = 0. This can be shown as follows. First, absence of radial acceleration D T u r = 0 
gives, upon using the conservation laws for E and J: 


1 - 


2 M 
~R 


2 - 


3 M 
~R 


mu 12 — ( 1 — 


3 M 
Jl 


mR V 2 = 2^(1- 


2 M 
~R 


E + 


JM 


R 


u\ (34) 


whilst the hamiltonian constraint (|32|) simplifies to 

■ t2 = 1 + R 2 u‘ p2 . 


2 M 


R 


u 


(35) 


These two equations can be solved for u t and u^ in terms of (R, E , J), implying that they 
are constant. An immediate consequence is, that S tr , £ r¥, and are constant as well, 
and actually vanishes. This follows directly from the absence of four-acceleration: 




du v 


M 


dr mR 3 

Then also the rate of change of Ti ttp must vanish: 


2 M 

7r” 


u 


■Tt* = 0 . 


(36) 


M , 

~R 1 


2 M\ dYP 


R 


dr 


, M\ 
^ |1 -« 


3 M 


R 


mdu* + 


JM 2 

~w 


u 


E 1 


2 M 
~R 


Now from eqs. (|34|) and (|35|) it follows that 

2 Af 


^u- 

m 


E = 2- 


3 M JM 


u 


* + R 2 u v2 . 


R J R mR 

These equations then allow the elimination of E and td, with the result that 


u* = 0 . 

(37) 

(38) 


JM (2M 


mR 2 \ R 


+ R 2 u v2 = Ru* 
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(39) 
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As for the total spin, for circular orbits the expression (33) can be written as 


/ 


-E ir 2 + 


R 2 E r ^ 2 

i 2 M 
1 R 


R 4 
M 2 



mu 


E 


+ 


(i-¥t 


[j-mifV] 2 . 


( 40 ) 


Thus for circular orbits u v and u f are constants which can be expressed in terms of R and 
J, in turn fixing E and / as well. 

6 Non-minimal hamiltonians 

So far we have studied the dynamics of compact spinning objects generated by the minimal 
geodesic hamiltonian H 0 . In this section we consider the non-minimal extension ([6]) 

K 

TT _ J? 

^ ^ j 


including the spin-spin interaction via space-time curvature. It is straightforward to derive 
the equations of motion: 


,-j.il _ 

=$■ 

7T/i = rng pu x u , 



TTfj = {7 T p , H} 

=> 

E R K \ p ~ 

^E kA E p(T V p R KXpa , 

(41) 

t iJU = {yt,h} 


D r TT = kE kA {Rj' a Z va 

- • 



Comparing again with the electro-magnetic force, the middle equation implies that in 
addition to the gravitational Lorentz force there is a gravitational Stern-Gerlach force, 
coupling spin to the gradient of the curvature. Therefore the coupling parameter k has 
been termed the gravimagnetic ratio [221113]. Like in the electromagnetic case ra the 
Pauli-Lubanski and Pirani-vectors are affected by this Stern-Gerlach force: 

D r S' J = | E pff V v R paaf) ) , 

(42) 

D r zV = -/cE kX RJ u Z v + (« + E^E kX R kXv(J u° - £- E^E^E p °S7 v R kXp(T . 


The second equation simplifies strongly for the special value 

1 

2 m 

In that case an initial condition = 0 is conserved up to terms of cubic order in spin. 










For the extended hamiltonian the conditions for the existence of constants of motion are 
modified. The total spin / defined in (13) is still conserved, but the conserved hamiltonian 
now is of course H = H 0 + H%. Finally we prove that the constants of motion J of the 
form (15) are preserved under this modification of the hamiltonian. To see this, observe 
that 

T 


( - a A V x R fwpa 


+ 


Vp(J 


(44) 


For the Killing-vector solutions (16) the right-hand side takes the form 


E^E P ° | _ a ^ xR/xvpa + fj llX R\ pa 


= - ^ u E pa (V^VpV* + VpV^V*) a. 


= 2 E '“' E '’"(V,.Vp + V„Vf)/ ? ™ = 0. 


(45) 


due to the anti-symmetry of the tensor f3 av . Therefore in particular the expressions (23) 
and (24) also define constants of motion in Schwarzschild space-time in the presence of 
Stern-Gerlach forces, as described by the non-minimal hamiltonian ([6j) . 

7 Conclusions 


In the context of general relativity the notion of point-masses is troublesome; any non¬ 
zero mass has a characteristic scale, typified by its Schwarzschild radius, describing its 
minimal size as defined by the corresponding horizon Therefore the approximation 
of a gravitating compact body as a point-like massive object in curved space-time requires 
the body to be small compared to the radius of curvature of the background space-time 
pa ee]. In addition, the mass must be small enough to ignore its effect on the space- 
time geometry at large. In the existing literature much effort has been put into obtaining 
effective equations of motion for compact objects by defining a position variable which 
can be interpreted effectively as the material center m Eg. One then computes the 
momentum and angular momentum in terms of a mass and momentum distribution in a 
finite neighborhood of this point. However, such a position variable is not unique, and 
moreover it often traces out a complicated world line, as shown for example by the well- 
known helical motion that is a solution of the Mathisson-Papapetrou-Dixon equations 

a edez!. 

In situations where a point-particle approximation of a spinning and gravitating body 
is appropriate, a complementary approach suggests itself by constructing a lagrangean or 
hamiltonian mechanics for a mass-point carrying spin in a curved space-time. In this letter 
we have chosen the hamiltonian point of view, as in our opinion this is most transparent 
in its results and application. In particular, the closed set of Dirac-Poisson brackets Q 
provides a unique and unambiguous starting point for the derivation of equations of motion 
for any representation of the spin degrees of freedom, allowing for a large class of physical 
implementations as fixed by the choice of hamiltonian. Two such choices, a minimal and 
a non-minimal one, have been presented and analyzed in this letter. 
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The minimal choice of hamiltonian is the one which also describes the geodesic motion 
of spinless particles. With this choice of hamiltonian the spin is covariantly constant along 
the world line, which is no longer geodesic due to spin-orbit coupling. It naturally provides 
a different implementation of the notion of position of the body, one for which now the 
Pirani vector Z is no longer taken to vanish. The advantage is that in terms of this choice 
of position variable the motion becomes tractable in non-trivial situations of practical 
interest; the motion in Schwarzschild space-time, as analyzed in sect. 5, provides a case in 
point. In addition, non-minimal hamiltonians can provide more complicated dynamics, as 
required for example for objects with non-vanishing gravimagnetic ratios [22] [33]. In this 
case the spin is subject to a kind of gravitational Larmor force, making it precess around 
field-lines of constant curvature. 

The question which effective hamiltonian to use for which physical system now becomes 
a matter of phenomenology. One should either derive the correct effective hamiltonian 
from first principles, connecting the formalism to the specific energy-momentum tensor, or 
determine it from experiments or observations. For the particular case of rotating black 
holes it could presumably be measured by observing gravitational waves from Extreme 
Mass Ratio binary systems involving a stellar-mass black hole; for a review see [48]. 
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